Spin kinetic theory 



— quantum kinetic theory in extended phase space 



Mattias Marklund, Jens Zamanian, and Gert Brodin 

Department of Physics, Umed University, SE-901 87 Umed, Sweden 

The concept of phase space distribution functions and their evolution is used in the case of en enlarged phase 
space. In particular, we include the intrinsic spin of particles and present a quantum kinetic evolution equation 
for a scalar quasi-distribution function. In contrast to the proper Wigner transformation technique, for which 
we expect the corresponding quasi-distribution function to be a complex matrix, we introduce a spin projection 
operator for the density matrix in order to obtain the aforementioned scalar quasi-distribution function. There 
is a close correspondence between this projection operator and the Husimi (or Q) function used extensively in 
quantum optics. Such a function is based on a Gaussian smearing of a Wigner function, giving a positive definite 
distribution function. Thus, our approach gives a Wigner-Husimi quasi-distribution function in extended phase 
space, for which the reduced distribution function on the Bloch sphere is strictly positive. We also discuss the 
gauge issue and the fluid moment hierarchy based on such a quantum kinetic theory. 

PACS numbers: 52.25.Dg, 51.60.+a, 71.I0.Ca 



I. INTRODUCTION 



In 1932, Eugene Wigner published the paper On the Quantum Correction For Thermodynamic Equilibrium yj] in which he 
set out to define the concept of a (quasi-)distribution function in order to approach the problem indicated in the title of his 
work. This seminal paper was very much the starting point for quantum kinetic theory (see also J2j-|4|]), in which the classical 
dynamical theory of Boltzmann was extended to the quantum domain. While the Wigner function is of interest in e.g. the 
interpretation of quantum mechanics I5l, l4(ill . the introduction of the quasi-distribution function (so called because of not being 
positive definite) of Wigner is not unique. In particular, within the field of quantum optics several other definitions naturally occur 
(see e.g. |6|] for an overview). These different definitions correspond to particular operator ordering and thus certain applications 
lends themselves naturally to different definitions. For example, when considering optical coherence normally ordered operators 
occurs naturally and hence the Glauber-Sudarshan P-distribution j4lll42ll is a convenient choice. On the other hand, anti-normal 
ordered operators i.e. the Q-function or the more general Husimi function l43ll are useful when dealing with quantum chaotic 
systems. For reviews of the subject see for example Refs. |44], HH] . Apart from the fields of plasma physics and quantum 
optics, where different phase space models have been applied with great success, the field of condensed matter physics and 
transport theory has also utilized the different kinetic(-like) approaches, such as semiclassical techniques [7], Green's function 
techniques @,|2,[n|,[l2], and diagrammatic techniques [10]. The research field known as quantum plasmas has recently become 
an area of intense investigations, with possible application to high intensity laser-plasma interactions [ 13, 20-22J1, hig h energy 
density physics [14], , nano- and sub micron-technology, e.g. quantum dots I23TI25I1 and plasmonic components |j26ll27ll . and 
nonlinear collective quantum problems lfl5l - [T9l I2H, l29ll . 

The above mentioned studies often assume an unmagnetized state of the system lf30h . However, significant influence can 
be exerted by external magnetic fields under the right conditions. For example, the magnetization dynamics of materials have 
important applications to data storage and memory production 13111 . Such dynamics is traditionally approached using the Landau- 
Lifshitz-Gilbert equation 13211 . Due to the large set of possible applications, as well as an interest to understand the basic 
dynamics of such systems, quantum plasma systems with spin has thus recently attracted a lot of interest from the research 
community (see e.g. IU9ll34ll and references therein). In the nonlinear regime spin solitons [36] and ferromagnetic behavior in 
plasmas can be found II37I1 . Up until this date, most of the presented studies have been of a theoretical nature. However, it is not 
difficult to envision future applications to e.g. plasmonic devices 12611 or femtosecond physics 13811 . 

In this paper we will give an overview of the field of spin extended phase space, and their implications for the formulation of 
a generalized kinetic theory. 



II. NONRELATIVISTIC MICROSCOPIC EQUATIONS: SCHRODINGER AND PAULI DYNAMICS 



A. The SchrSdinger description 



The dynamics of an nonrelativistic scalar electron, represented by its wave function if/, in an external electromagnetic potential 
' is governed by the Schrodinger equation 

d\li h 9 

ih^ + —V 2 ifr + e<f>i/t = 0, (1) 
at Lm e 
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where fr is Planck's constant, m e is the electron mass, and e is the magnitude of the electron charge. Using the decomposition 
ip = yfn expiS /ft into the amplitude and phase IB9I1 . (1) can be written as two real equations. Here n is the amplitude and S the 
phase of the wave function, respectively. Thus, Eq. (HJ becomes 

^ + V ■ (nv) = 0, (2) 
ot 

and (after taking the derivative of the equation for the phase) 

m e \— + v- V v = eV0 + - — V — ^- I , (3) 



dt ) 2m e \ ^Jn 

where the velocity is defined by v = VS /m e . The last term of Eq. © is the gradient of the Bohm-de Broglie potential, and is 
due to the effect of wave function spreading, giving rise to a dispersive-like term. We also note the striking resemblance of Eqs. 
© and © to the classical fluid equations. 



B. The Pauli description 

Through the Dirac Hamiltonian 

H = ca ■ (p + eA) - e$ + /3m e c 2 , (4) 

the spin is introduced. Here a - (ai,a2, 0:3), e is the magnitude of the electron charge, c is the speed of light, A is the vector 
potential, <p is the electrostatic potential, and the relevant matrices are given by 

<x \ „ / I 



rrOV '= 0-1 (5) 



Here I is the unit 2x2 matrix and cr = (o"i , 0"2, (T3), where we have the Pauli spin matrices 



o)' cr2 = (° j' and ^3=(J \\ (6) 
From the Hamiltonian the weakly relativistic limit (in v/c) gives 

H = t— (p + eA) 2 + ^-B cr — ecp. (7) 

We see that the electron's magnetic moment is given by m — -yLB{ii\<^W)l{^W), where fig - eh/2m e is the Bohr magneton, 
giving a contribution -B • m to the energy. The latter shows the paramagnetic property of the electron, where the spin vector 
is anti-parallel to the magnetic field in order to minimize the energy of the magnetized system. Given any operator F and the 
hamiltonian (©, the relation dF/dt = dF/dt + (l/ih)[F, H], where [, ] is the Poisson bracket, gives the evolution equations for 
the position and momentum in the Heisenberg picture [ 39] 

$ = — (p + eA) = v, (8) 

dt m e 



m e — = -e (E + v x B) - -fx B V(B ■ s), (9) 

dt n 

respectively, while the spin operator satisfies the evolution equation 

ds 2 

- = -^8X8, (10) 
dt n 

showing the spin precession in an external magnetic field. Here the spin operator is given by s = {fi/2)cr. The above equations 
thus gives the quantum operator equivalents of the equations of motion for a classical particle, including the evolution of the spin 
in a magnetic field. 
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The non-relativistic evolution of spin-^ particles, as described by the two-component spinor !f( a ), is given by the Pauli equation 
(see e.g. JH) 



H I ie . \ 2 
— V + -A -, B B-cr + , 



<A = 0, (11) 



where A is the vector potential, fig — eh/2m e is the Bohr magneton, and <x = (cr\ , <x 2 , 1x3) is the Pauli spin vector. 

We may, as in the Schrodinger case, decompose the electron wave function iff into its amplitude and phase. Thus, we let 
iff = yn exp(/5 /K)<p, where tp is a unit two-spinor. Then we have the conservation equations 

— + V-(nv) = (12) 
ot 



and 



m e ^- = -e(E + v x B) + (^^| - ^(V ® B) ■ s - —V • (nZ) (13) 



respectively. The spin contribution to Eq. (fT3l is consistent with the results of Ref. |70!]. Here the velocity is defined by 

1 / \ 

v = — I VS - iHip^Vip) + , 

the spin density vector is 



(14) 



h + 

s = -ip'o-ip, (15) 

which is normalized according to |s| = h/2, and we have defined the symmetric gradient spin tensor 

Z = (ViJ®(V/). (16) 
Moreover, contracting Eq. (fTTT i by ificr, we obtain the spin evolution equation 

ds _ {2(i B _ B _ J_ [5a(n5 a s)] | x s (17) 



dt ( H m e n 

We note that the last equation allows for the introduction of an effective magnetic field B e ff = (2/i#/ft)B - (m e ny x [d n (nd a s)]. 
However, this will not pursued further here (for a discussion, see Ref. l39lo . 

III. COLLECTIVE PLASMA DYNAMICS 

As pointed out in the previous section, the route from single wavefunction dynamics to collective effects introduces a new 
complexity into the system. At the classical level, the ordinary pressure is such an effect. In the quantum case, a similar term, 
based on the thermal distribution of spins, will be introduced. 

A. Fluid model 

1. Plasmas based on the Schrodinger model 

Suppose that we have electron wavefunctions, and that the total system wave function can be described by the factorization 
iff(x\,X2, . . .Xjv) = • ■ • i'N- For each wave function \ff a , we have a corresponding probability V a . From this, we first define 
iff a = n„ expiiSn/K) and follow the steps leading to Eqs. (O and (01. We now have such equations the wave functions {iff a }- 
Defining H 

n=yp a n a (18) 
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and 

N 



v = (\ a ) = > J , (19) 



we can define the deviation from the mean flow according to 

w ff = Vqt — V. (20) 
Taking the average, as defined by ( fT9l ), of Eqs. (O and ([3]l and using the above variables, we obtain the quantum fluid equation 

J+V-(mO = (21) 

and 

^(| + v.v)v = e „V,-Vp + £v((^)), (22) 

where we have assumed that the average produces an isotropic pressure p = m e n(\v/ a \ 2 ) We note that the above equations still 
contain an explicit sum over the electron wave functions. For typical scale lengths larger than the Fermi wavelength Af, we may 
approximate the last term by the Bohm-de Broglie potential [ 80] 

^\..v»vs (23) 



/ V» 

Using a classical or quantum model for the pressure term, we finally have a quantum fluid system of equations. For a self- 
consistent potential <p we furthermore have 

V 2 = -(" - «;)■ (24) 



2. Spin plasmas 

The collective dynamics of electrons with spin, based on a fluid model, was presented in Ref. 1 34] . For the present discussion, 
we will follow Refs. I34ll35ll . Suppose that we have Af wave functions for the electrons with magnetic moment -ps, and that, as 
in the case of the Schrodinger description, the total system wave function can be described by the factorization i// = i//\i//2 . . ■ ^n- 
Then the density is defined as in Eq. ( fTBl and the average fluid velocity defined by ( fT9b . However, we now have one further fluid 
variable, the spin vector, and accordingly we let S = (Sq.). From this we can define the microscopic microscopic spin density 
S a = s a - S, such that (S a ) = 0. 

Taking the ensemble average of Eqs. ( fT2b we obtain the continuity equation d2"Tl i, while the ensemble average applied to ( fT3b 
yield 

id \ h 2 n (V 2 tJH\ 

meU [jt + V ' V j V = ~ m (E + V X B) " VP + 2n7 V [~^hTj + Fspi " (25) 



and the average of Eq. ( [T7| > gives 



dt / h 



n\— +v-V|S= -^-BxS- V-K + n spin (26) 



respectively. Here the force density due to the electron spin is 

F spin = -3^(V«B) ■ S - — V ■ Ul + 1)1 - —V • [MVS a ) ® (VS a a ) + n{VS a a ) 9 (VS a )], (27) 

consistent with the results in Ref. lf70ll . while the asymmetric thermal-spin coupling is 

K = n(w a ® S a ) (28) 
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and the nonlinear spin fluid correction is 

Q S pi n = — S x [d a (nd a S)] + — S x [daini&'Sa))] 
m e m e 

+ILl^L x {d a [n a d a (S + S a )]}) (29) 
m e \n a J 

where E = (V5 fl ) ® (VS fl ) is the nonlinear spin correction to the classical momentum equation, E = ((VS(a)a) ® C^L))) i s a 
pressure like spin term (which may be decomposed into trace-free part and trace), and [(V ® B) • S ]" = {d a Bb)S b . Here the 
indices a,b, . . . = 1, 2, 3 denotes the Cartesian components of the corresponding tensor. We note that, apart from the additional 
spin density evolution equation d26l ). the momentum conservation equation (l25b is considerably more complicated compared to 
the Schrodinger case represented by d22l) . Moreover, Eqs. ( f25l ) and ( f26b still contains the explicit sum over the N states, and has 
to be approximated using insights from quantum kinetic theory or some effective theory. 

The coupling between the quantum plasma species is mediated by the electromagnetic field. By definition, we let H = 
B/yizo _ M where M = -2n/i#S/ft is the magnetization due to the spin sources. Ampere's law V x H = j + £od,E takes the form 

VxB^„G+VxM) + iJ, (30) 
c z at 

where j is the free current contribution The system is closed by Faraday's law 

VxE = — . (31) 
at 



IV. DISTRIBUTION FUNCTIONS IN QUANTUM MECHANICS 

The Wigner distribution 11|] function is the quantum version of the classical phase-space distribution function. It was developed 
during the 20s and 30s following the success of classical non-equilibrium statistical mechanics. However, since in the quantum 
case the Heisenberg uncertianty principle has to hold, the introduction of a distribution function is not straightforward. There 
are in fact infinitely many different ways to construct a quantum analogue to the classical distribution function. Except for 
the Wigner function other popular distribution functions are the Husimi (or Q-) function B43I1 and the Glauber-Sudarshan P- 
distribution EI HI. 

There are some properties that are natural to impose on the distribution function: 

• Non-negativity, i.e. /(x, p) > 0. 

• Correct marginal distributions, i.e. 

/(x) = J d 3 pf(x,p) (32) 

/(p) = J d 3 xf(x,p), (33) 

should yield, respectively, the probability density to find a particle around the space point x or with a momentum around 
the point p in momentum space. 

Other properties that can be imposed on the distribution function are that it is real valued, linear in the density matrix and that 
one can find a complete orthonormal set in the space the distribution function. Since the particle position and momentum cannot 
be simultaneously known, it is not possible to find a quantum distribution function satisfying all of the conditions above. In order 
for the description to be complete and equivalent to the density matrix formalism it must be possible to calculate the expectation 
value of any observable using the distribution function. To do so the operator corresponding to the observable is first mapped 
to a phase space function O = 0(x, p) — > 0{x, p) and the phase space average weighted by the distribution function is them 
calculated according to 



«s>=r„wc*.P>«,P>. 04) 

The mapping from the operator space to phase space depends on which distribution function is used, see ll45tl for details. In this 
proceedings the phase space distribution function we will use is the Wigner distribution function OJJ] . Given the density matrix p 
it is defined by 

f w (x, p) = — i-j f d 3 y e- ip » n p(x + y/2, x - y/2). (35) 
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It produces the correct marginal distributions, but it may become negative in some regions. The correspondence between opera- 
tors and phase space functions is in this case is called Weyl correspondence ||2|] and is given by 



0(% p) <^ 0(x 



p) = J d'ye 



-i'p-y/fi 



<x + y/2|0|x-y/2>. 



(36) 



The discussion so far has been in terms of phase space distribution functions, however, there are similar ways to construct 
distribution functions for particles with spin E 56T - D9I1 . Here the spin analog to the Q-function will be utilized. It is defined by 



/(§) = Tr 



— (1+g-o-) p 

4-7T 



(37) 



where s is a unit vector which will have the same role as the phase space variables x and p, cr denotes a vector with the Pauli 
matrices [Eq. ©J as its components and p is the 2x2 density matrix. Similarly to the phase space case the spin distribution can 
be used to calculate the expectation value of any operator acting on the spin degree of freedom by integration. Since the Pauli 
matrices satsifies cr? = /, for i = 1, 2, 3 where / is the 2x2 unit vector the only two operators that can act on spin degree of 
freedom are the unit operator / and the Pauli matrices cr. The expectation values of these are given by 



</> = j d 2 sf(s) 
(cr) = 3 j d 2 ssf(s). 



(38) 
(39) 



The somewhat peculiar occurrence of a factor 3 can be understood by noting that a pure quantum state is "smeared out" over the 
whole sphere in contrast to the classical case where a dipole moment may point in a given direction, see l29Tl for details. 



A. Evolution equation 

The time evolution of the density matrix is given by 

mf = [&,p}. 



(40) 



For non-relativistic electrons the Hamiltonian is given by Eq. (0. Taking this as the starting point and using the two transforma- 
tions Eqs. d35l l and ( l37b the evolution equation for the spin distribution function is obtained 



dt J 



— (E + v x B) + — V V [(V S + s) • B) 



< V v / - ^(§ x B) • V s f 
n 



^ ( y_ v . A) _^( B .^ + s.B)lfcsin(fv,.V l ,)-V,.V l , 

m e m \ / J [ n \ 2m j 



f 



(41) 



-A • V, + —[(A • V X )A] • V v - -^(§ x B) • V, 

171 YYl n 



2m e 



V,- V,, -1 



/■ 



where the Coulomb gauge and the variable transformation m e \ = p + eA(x) has been used to simplify the equation somewhat. 
In the equation above the function of the operators are defined by their Taylor expansions. A more appealing version of the 
evolution is obtained if the long scale length limit is considered. Neglecting terms of order h 2 and higher the evolution equation 
is written as 



dt J 



— (E+vxB) + — V t (s-B) 



V v / - ^(§ x B) • Vj/ - —[V X (B ■ V,)] • V v / = 0. 
n m e 



(42) 



Note that the last term contains derivatives both with respect to the velocity v and the spin s. It is interesting to note that we 
can derive a similar equation on semiclassical grounds. Starting with the distribution function F(x, v, s) and setting the total 
derivative along particle paths to zero, we obtain 



„ DF dF „ d\ „ ds „ 

= = — + v ■ V X F + — ■ V V F + — • V SJ F. 

Dt dt dt dt 



Now using equations (9) and (10) we obtain 



dt XJ 



— (E + vxB)+ — V t (s-B) 

m e m e 



V v /-^(§xB)-V s / = 0. 
n 



(43) 



(44) 
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This equation has already been studied in ll64ll . It is there shown to give rise to new oscillation modes due to the anomalous 
magnetic moment of the electron. A similar equation has also been studied in lf75tl where it is investigated whether spin may 
be of importance in magnetic confined fusion experiments. The difference between Eq. d42l and the semi-classical case Eq. 
(l44l is due to the fact that the quantum mechanical probability distribution is smeared out compared to the classical distribution 
function. 



B. Charge and current density 

To obtain a self-consistent mean-field description of a plasma the evolution equation for the distribution function has to be 
coupled to Maxwell's equations. This is done by noting that the quantum mechanical operator corresponding to the density is 
«(x) = <5(x - x) which according to Eq. (l36i l is transformed into the phase space function n(x) = 6(x - x'). Hence the electron 
charge density is given by 



-en M = -e J M^sSix - x')f(x' v,M) = -e j>v/(x, v, s, 0. 



(45) 



To this one should add the corresponding charge density for the positive species or add a neutralizing homogenous stationary 
background charge eriQ, where no — J cPxcPvd 2 sf . Due to the spin magnetic moment there will be a magnetization current 
density and the total electron current density is given by |p70h 

j e (x,t) = -e J d 3 vd 2 svf(x,v,s, t) - 3// B V. v x J d 3 vd 2 ssf(x,\,s, t). (46) 
The first term is seen to be the the free current and last term the curl of the magnetization due to the electron spin. 



V. IMPLICATIONS OF GAUGE INVARIANCE 

The definition of the Wigner function d35l l is not gauge invariant since it is a function of the gauge dependent canonical 
momentum rather than the gauge independent kinetic momentum m e v — p + eA(x). The theory above is hence only valid in the 
Coulomb gauge. It is possible to modify the definition to obtain a gauge independent Wigner function [68]. In principle, there is 
nothing that prevents us to use a gauge dependent Wigner function as long as care is taken when doing gauge transformations. 
However, problems may arise when calculating for example the second order moment of the velocity (vj$/). One might then be 
tempted to calculate 

J d 3 x d 3 WiVj f(x, v, t) = J d 3 x d 3 p [pi + eA,(x)] [pj + eAj(x)] f(x, p, ?)■ (47) 

However, the phase space function which is related to the operator v,Vj is not [/?, + eAi(x)][pi + eAj(x)]. In order to obtain the 
right function it is necessary to first put operator [pj + eAj(x)][pj + <?A ; (x)] in Weyl-ordering [2] and then make the substitution 
x — > x, p — > p. This is in general difficult to do since the vector potential is a function of x. In the current proceedings we are 
only considering the charge and current densities Eqs. d45b and ( l46b and this problem never arises. 



A. A gauge invariant distribution function for the case of spin-extended phase space 



For completeness the fully gauge invariant distribution function is given here. Following Ref. 16811 the Wigner matrix is 
defined by 



W GI (x 



' v ' a - fi - f) = (2W I dhexp H? v ' [ z ' e Li dTA(x + TZ ' °]} p ( x + Z r a '' x ' l' fi ) • 



(48) 



where v is the velocity. The explicit dependence of the vector potential in this construction is there to compensate for the phase 
factor which the wave function acquires under a gauge transformation. Using the spin transformation ( l37l i. we obtain a fully 
gauge invariant distribution function. The evolution equation for the gauge invariant distribution function without spin was 
derived in Ref. B68I1 . It is straightforward to generalize this equation to include spin with the result 



rGI 



" r + (v + Av) ■ V,/ GI — — f(v + Av) x B + £l • V r / GI - ^ V,[(8 + V,) • B] • V,,/ 01 - ±Z \s x (S + AS)1 • V. S / GI = 0, (49) 



dt 
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where we have defined 

-1/2 



- f 1/2 / ihT \ 
E= drE x+ — V v =E( 

J-i/2 \ w e ; 

- r 1/2 / i»r \ 

B= drB x + — V v =B(x) + b( 

J- 1/2 \ w e ; 



x) + e(x) 
x) 



1/2 



Av 



iqh f 

m 2 J_i/2 
ft f 1/2 



<irrB x + — V,, X V,, = 7 



B(x) c/TTsin — V c ■ V„ 

J-l/2 W / 



/-1/2 
-1/2 



X Vv= - 



eft 2 
12mf 

h 2 



(50) 
(51) 

B(x)x V v (V v ■ V v ) + 0(h 4 i)52) 



i% C' I ifa \ <- -» ft f ' /-rft ^ ^ ft 2 <- -» 

AB = </ttB x + — V v V, • V v = — B(x) rfTTsin — Vv ■ vJv, • V v = — -rB(x)(V* • V,) 2 + 0(h 4 ), (53) 

m c J-i/2 \ '« / J-i/2 \m e I 12m| 



and 



r l/2 

E = E(x) rfr 

J-l/2 
-1/2 



(tH 




V„J - 1 


cos 


V,- 


\ m 






(tTi 




V v j - 1 


cos 


V. r 


\ to 







tl 2 



24m 2 

h 2 

'24m 2 



E(x)(V,- V v f+0(h 4 ) 
B(x)(V,- V,) 2 +CP 4 ), 



(54) 
(55) 



so that (cf. (40)) 
d/ GI 



+ v ■ V,/ GI - — Ie + vxB - ^ V a -[(s + V.0 ■ B] 1 ■ V,,/ GI -^(sxB)- V. ; / GI 
m e { m e J ft 



Av • V.J GI + — + (v + Av) x b + Av x B - — V x [(s + V. f ) b] } ■ V v / 



,GI 2 Vb 



[s x (b + AB)] • V. ; / GI , 



(56) 



or, to lowest order in ft, 
df 



^ + v ■ V,/ U1 --{E + vxB- ^V,[(s + Vj) B] \ • V,,/ GI - (§ x B) • V,/ 
of to 



2he 



24m 2 



E + vxB- ^V. r [(§ + V. ; ) • B] ■ V v (§XB) • V ; (V.v • Vv) 



ft 



(57) 



2e 



2e 2 



+ — (Bx V„> V, T 



(Bx V„) x B 



• Vv 



(V. t ■ Vv)/ 01 , 



The gauge invariant Wigner function has a modified Weyl correspondence which is well suited for calculating fluid moments. 
In order to obtain the phase space 0(x, v) function which corresponds to an operator 0(x, \), all products of the operators x and 
v = [p + eA(x)] /m e are first ordered in a symmetric form using the commutation relation x and p and then the substitution x — » x 
and v — > v is taken (details can be found in jUt ). 

The transformation between the gauge-invariant extended Wigner function and the gauge dependent dito can be obtained 
through a kernel function or through an differential operator defined through 



/ GI (x, v, t) = h(x, -i(fi/m e )V v , i)f(x, v, t), 



where 

ft(x, u, t) = exp[/6>(x, u, f)] = exp 



(-a M+ £;, 



— u-|-A(x, t)+ I dT A(r,x + tu)| 





ie 1 2 sinh z 


- l)A(x,0 


| = exp 


[~n u i z 




c=uV> / 



(58) 



(59) 



= 1 + [exp(/0) - 1] = 1 + ft(x, u, t), 



where the part ft gives the QM correction to between the gauge-invariant and gauge-dependent distribution functions. Thus, 
expressed through the gauge dependent extended Wigner function, Eq. d56*l l contains terms proportional to h 2 also on the left 
hand side. Therefore, what is meant with quantum mechanical corrections is dependent on the corresponding choice of Wigner 
function. We explicitly have 



/ GI (x, v, t) = exp 



A(x, t) 

m„ 



2 sin((ft/TO e ) V* ■ Vv) 



- 1 



Vv 



eft 2 



(h/m e ) V. t • Vv 

giving us the relation between the gauge independent and gauge dependent distribution functions as an infinite series. 



/(x, v, t), = /(x, v, f) - T f T A(x)- V,. (V* • Vv) 2 /(x, v, t) +0(h 4 ), 

(60) 
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VI. SUMMARY AND DISCUSSION 

In the present paper we have given an overview of the field of quantum kinetic evolution equations, as given by Eqs. ( HTt 
and d56l l, for a quasi distribution function of electrons, based on a Wigner transformation of the density matrix, together with 
a spin operator contracting the 2x2 Wigner-matrix to a scalar function /(x, p, s) (or, equivalently, f(x, v, s)). The free current 
and the magnetization can be directly computed from the quasi distribution function, and hence Eq. fiTI ) (or the gauge invariant 
alternative, Eq. j49l ) together with Maxwell's equations with the sources d45l l, and d46b , form a closed set. A discussion of the 
gauge problem was given. 

The gauge independent distribution can be used to derive fluid moments i8ll l82tl . The resulting equations can be used to 
investigate for example nonlinear phenomena. However, even in the linear case it is necessary to keep up to the sixth order 
moment in order to retain the lowest order quantum contribution. Hence, at least for linear analysis the kinetic equation is a 
better tool. 
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